Nonadiabatic Superconductivity and Vertex Corrections in Uncorrelated Systems 



O 

o 

(N 

-t— > 
o 

o 

o 

m 

o 
o 

u 

Oh' 



X3 

o 



o 



I 

o 
o 



X 
S3 



P. Paci^, E. Cappelluti^, C. Grimaldi-^ and L. Pictronero^ 
^ Dipartimento di Fisica "A. Volta", Universitd di Pavia, Via Bassi 6, 27100 Pavia, Italy 
and Istituto Nazionale Fisica della Materia, Unitd di Roma 1, Italy 
^ Dipartimento di Fisica, Universitd di Roma "La Sapienza" , Piazzale A. Moro, 2, 00185 Roma, Italy 
and Istituto Nazionale Fisica della Materia, Unitd di Roma 1, Italy 
Ecole Polytechmque Federale de Lausanne, Departement de Microtechnique IPM, CH-1015 Lausanne, Switzerland 

We investigate the issue of the nonadiabatic superconductivity in uncorrelated systems. A local 
approximation is employed coherently with the weak dependence on the involved momenta. Our 
results show that nonadiabatic vertex corrections are never negligible, but lead to a strong sup- 
pression of Tc with respect to the conventional theory. This feature is understood in terms of the 
momentum-frequency dependence of the vertex function. In contrast to strongly correlated systems, 
where the small q-selection probes the positive part of vertex function, vertex corrections in uncor- 
related systems are essentially negative resulting in an effective reduction of the superconducting 
pairing. Our analysis shows that vertex corrections in nonadiabatic regime can be never disregarded 
independently of the degree of electronic correlation in the system. 
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One of the characteristics of the high transition tem- 
perature supercoriductors (HTSC) is the small density 
of charge carriers As a consequence, the scale of the 
electronic dynamics, the Fermi energy Ep, is compa- 
rable with the typical phonon frequencies Wph and the 
adiabatic parameter Upi^/Ep becomes significant. This 
situation opens the way to a scenario in which nonadi- 
abatic effects are relevant. From a diagrammatic point 
of view, nonadiabatic effects can be taken into account 
by the inclusion of vertex corrections arising from the 
breakdown of Migdal's theoremE However, in general, 
this amounts to consider an infinite set of nonadiabatic 
diagrams whose resummation is a formidable task. Nev- 
ertheless, it is possible to formulate a perturbative theory 
by assuming that the ordej of magnitude of the first ver- 
tex corrections, Xcjpi^/ Epp is-small enough to be treated 
as an expansion parameter .Ba In principle, this assump- 
tion is fulfilled by nonadiabatic weak coupling systems 
(A < 1 and uiph/Ep 1) or moderately nonadiabatic 
strong coupHng materials (A ~ 1 and uiph/Ep < 1). 

In previous studies we have shown how such a per- 
turbative approach accounts for some of the anomalous 
properties of the HTSC-materials. In particular, large 
values of Tc, compared with the ones predicted by the 
Migdal-Eliashberg (ME) theory, are related in a natu- 
ral way to the opening of nonadiabatic channels in the 
Cooper pairing with no need of assuming large values of 
A.LfB To understand from a microscopical point of view in 
which way the nonadiabatic channels affect the supercon- 
ducting properties, and in particular the critical temper- 
ature Tc, a detailed study of the momentum- frequency 
structure of the vertex function is needed. In fact, ver- 
tex function presents a complex behaviour with respect 
to the momentrmi, q and the frequencies u of the ex- 
changed phonon.Bcl In particular the vertex function has 



been shown to be positive for small values of q and neg- 
ative for large values of q compared to w, leading re- 
spectively to an enhancement and to a decrease of the 
superconducting pairing. The evaluation a priori of the 
nonadiabatic effects on Tc is thus not at all easy, since it 
will depend on the total balance of negative and positive 
parts of the vertex function. In this perspective, specific 
properties of real materials become very important, since 
they can modify the balance favouring or disfavouring 
positive or negative parts and determining an enhance- 
ment or a suppression of the critical temperature. 

In particular, in strongly correlated systems the elec- 
tron correlation due to onsite Coulomb repulsion has 
shown to be actually responsible for a predominance of 
small-q scattering yielding-M effective modulation of the 
electron-phonon coupling.QEi In this situation the posi- 
tive part of the nonadiabatic vertex corrections is mainly 
probed, leading to an net increase of the coupling in 
Cooper channel and to a corresponding enhancement of 
the critical temperature Tc.oQ On the other hand in un- 
correlated systems the q-dependence is weak and the neg- 
ative part of the vertex corrections at large momenta is 
expected to lead to a resulting decrease of the Cooper 
pairing and of Tc. 

The aim of this short communication is twofold. First, 
we show that nonadiabatic effects cannot be neglected 
also in uncorrelated materials (structureless electron- 
phonon interaction) and vertex corrections play a pri- 
mary role in suppressing the superconducting pairing as 
long as ujph/ Ep is not negligible. Second, we argue that 
the nonadiabatic superconductivity developed by us for 
the small q scattering regime is also a rather good ap- 
proximation of the uncorrelated case, providing therefore 
a unified and reasonable description of both the corre- 
lated and uncorrelated cases. 
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In coQventional metals, according to Migdal's 
theoremjj the smallness of the adiabatic parameter 
uiph/Ep permits to describe successfully the electron- 
phonon coupled system by neglecting the vertex cor- 
rections in the electron-phonon interaction. The ap- 
plication of Migdal's theorem to the superconducting 
state has led to the ME equations of superconductivity,Ll 
which accuratelyi-dfiscribe the properties of conventional 
superconductors.t2l A different situation is encountered 
in nonadiabatic materials, as we have briefly discussed 
above, where the breakdown of Migdal's theorem is ex- 
pected. The relevance of the nonadiabatic corrections 
can be established by evaluating the vertex function 
schematized as 

P(w„, LUra; q) = 5^ X! ~ ^l)G{p, LUl) 

xG{p + (i,uji+uj„i) , (1) 

where D and G represent respectively the phonon and 
the electron propagators, g is the electron-phonon matrix 
and LUi are fermionic Matsubara 



0J„ 



element, and w„ 
frequencies. We adopt a simple Einstein spectrum with 
frequency ujq for the phonon propagator D{u}n — lui) = 

As shown by Migdal in his pioneering workp the vertex 
function, given in the Eq. (|l|), scales as XiOph/Ep where 
A = 2Nog'^/ujQ is the dimensionless electron-phonon cou- 
pling and A'o is the electronic density of states (DOS) 
at the Fermi level. To obtain this result typical phonon 
momenta are assumed to be of the order of the Debye mo- 
mentum ^ kp. However, as briefly above discussed, 
this assumption breaks down in strongly correlated sys- 
tems where the predominance of small-q scattering is 
important to establish the enhancement of the critical 
temperature. 

The small-q selection due to strong Coulomb repul- 
sion in strongly correlated systems can be simulated by 
a cut-off qc in the exchanged momentum space so that 
^ X]q ^('Zc — |q|) which restricts the momentum inte- 
grations. In this situation we replace the vertex function 
P{ijJn, ^m] q) by its average over momenta P{LUn, ujm', 9c) 
that depends only on the frequencies and on the cut-off 
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We can generalize the ME equations to include the mo- 
mentum average of the first order vertex corrections due 
to the breakdown of Migdal's theorem. The generalized 
ME equations in the nonadiabatic regime for the self- 
energy renormalization function ^ and for the supercon- 
ducting gap A are the following:! 



Z{uJn) 



1 



LJn, OJr, 
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where rjm = 2 a,Tcta,n{Ep /[Z (cjm) ^m]} and Qc = qc/2kp 
is the dimensionless cut-off. The kernels of the equa- 
tions (^) and are respectively given by: 

r2(w„, Wm,Qc) = AI?(W„ - Wm)[l + XP{u)n,U)m,Qc)\, 



V^{uJn,UJm, Qc) = XD{uJn " W„)[l + 2AP(w„ , W„ , Qc)] 
+ X'^C{uJn,UJm,Qc)- 

Explicit expressions of the vertex P and cross C func- 
tions have been obtained analytically for small values of 
Qc in Refs. ^,|ll]. Higher order nonadiabatic corrections 
in Tz and Fa, here not taken into account, should be 
explicitely included in th£_|extreme nonadiabatic regime 
where XP ~ 1 (AC ~ 1).E3 This intriguing issue is how- 
ever beyond the purpose of our paper, and all the follow- 
ing results apply only when Xujph/Ep is small enough to 
permit truncation of higher order vertex corrections. 

The nonadiabatic equations of superconductivity 
Eqs. (H) and (^) have been numerically solved and the 
resulting Tc has found to be strongly enhanced with re- 
spect to the adiabatic case due to the inclusion of vertex 
corrections and to the presence-of electronic correlation 
that favours small-q scattering.Bi3 

A different situation is encountered when we look at 
uncorrelated systems where the q-dependence of the rel- 
evant physical quantities is weak. In this case no re- 
striction in the momenta integration is expected and the 
momentum average of the vertex function becomes an al- 
most exact approximation. As easily seen from Eq. (|^) 
this corresponds to a local theory in which the vertex 
function becomes 

Pioc(Wn, W,„) ^ PiuJn, ^Jm', q) 

q 

= 5^ X! ~ ^l)Gloc{^^l)Gloc{l^l + LUm), (5) 

where Gioc is the local electron propagator: 
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de- 



N{e) 



'iuJnZ{uJn) 



For a direct comparison we adopt the same simplifica- 
tions of Refs. namely we considered a half- filled con- 
stant DOS band with N{e) = Nq, -Ep < e < Ep. Ep 
represents then the Fermi energy. 

The nonadiabatic equations of superconductivity for 
uncorrelated systems are thus formally obtained by sub- 
stituting the local vertex function given by the Eq. (H) 
in the kernels Tz and Fa of the equations (||) and (|j). 
The numerical solution of such equations in local regime 
follows the usual scheme. Therefore, without entering in 
details, we are going to discuss the results. 

In Fig. I^a we show the superconducting transition 
temperature of uncorrelated nonadiabatic systems as 
function of dimensionless electron-phonon coupling A for 
three different adiabatic parameters. To evidence the 
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FIG. 1. (a) Tc as function of A in uncorrelated nonadi- 
abatic systems. Solid lines correspond to the non crossing 
approximation; dashed lines correspond to the nonadiabatic 
local theory. Both the cases are shows (from top to bottom) 
for ujq/Ef ~ 0.1, 0.4, 0.7. (b) Behaviour of effective pairing 
interaction VbB as function of A in non crossing approxima- 
tion (solid line) and in the nonadiabatic local theory (dashed 
lines) . 

effects of the vertex corrections, the results of the nona- 
diabatic local theory are compared with those obtained 
by the conventional one for values of A up to the un- 
physically large A = 3, where of course the perturbative 
approach breaks down {\ujq/Ep ~ 1). Solid lines repre- 
sent the ME solutions where nonadiabatic effects are just 
included by considering finite energy bandwidth. This is 
thus equivalent to a non crossing approximation (NCA). 
Dashed lines are the numerical results of Eqs. (||)-(|^), 
where the vertex function is given in the local theory by 
Eq. (H). From top to the bottom, different lines corre- 
spond to adiabatic parameters: LOajEp — 0.1, 0.4, 0.7. 

As shown in Fig. |l|a nonadiabatic effects in uncorre- 
lated systems lead to a drastic reduction of the critical 
temperature with respect to the conventional theory. 
This result confirms the above qualitative discussion sug- 
gesting that in uncorrelated systems the negative part of 
the vertex function is more relevant yielding therefore an 
effective reduction of the pairing. The increase of this 
effect by increasing both the electron-phonon coupling A 
and the adiabatic parameter loq/Ef seems also to point 
towards a similar conclusion, since the vertex corrections 
scale as Xcuo/Ep. 

In order to quantify this concept we parametrize the 
magnitude of the "effective" superconducting pairing, 
Vcs, with the static limit of the superconductive ker- 
nel: Vcs — rA(t^n = 0,uj,n = 0). In Conventional ME 
theory Vcs reduces to the simply bare electron-phonon 



coupling constant VcS = A, while the opening of nona- 
diabatic channels strongly affects Ves through the ver- 
tex function calculated in the local theory. The depen- 
dence of the effective interaction Vcs on the bare electron- 
phonon coupling A in both the cases is plotted in Fig. 
The similarity between the behaviours of Tc and of Vcs 
is striking pointing out that this particular limit of the 
superconductive kernel and of the vertex function is di- 
rectly reflected. on the critical temperature T^. the static 
negative limitEl of the vertex function reduces Vcg and 
induces a strong supression of Tc- The non monotonic 
behaviour of Tc in uncorrelated nonadiabatic systems is 
thus simply related to the non monotonic underlying be- 
haviour of VeS'- the maximum is determined roughly by 
the value of A which gives dVcS /dX = corresponding to 
4AP/oc(i^n = 0,a;m = 0) = — 1 (cross function is negligi- 
ble and can be omitted for the discussion). We conclude 
that nonadiabatic effects are dominant also in uncorre- 
lated systems where they lead to a strong reduction of 
Tc. 

Let us address now whether the nonadiabatic the- 
ory developed in Ref. for small momentum transfers 
(small Qc) are capable of reproducing the results of the 
local theory {Qc = 1). In fact, it would be interesting to 
define a common approach which would permit to span 
from uncorrelated to correlated cases. From this point 
of view, the small-Qc expansion of the vertex function 
used in Refs. appears a promising tool to extrapolate 
to intermediate QcS. In particular, as we can see from 
Eq. (H), the uncorrelated case corresponds in this frame- 
work to a Qc = 1- We stress that the equivalence between 
this procedure and the local theory is only formal since 
using Qc = 1 in a small-Qc expansion is obviously an 
approximation. 

In Fig. ^ the critical temperature Tc and the "effec- 
tive" superconducting pairing VeS calculated within the 
small-q approach with Qc = 1 are shown. The overall 
features of both Tc and VcS are quite the same as in the 
local theory with a slight overestimation of Tc and Vcg. 
This is not unexpected since the small-q expansion em- 
phasizes the positive region of the vertex function and 
as a consequence underestimates the net negative mag- 
nitude of it. However the qualitative agreement between 
Figs. |l] and || is quite good suggesting that the nona- 
diabatic evaluation of the vertex function based on the 
small-q expansion can actually interpolate from weak to 
large correlation. This is confirmed in Fig. |^ where we 
plot the evolution of the critical temperature Tc in nona- 
diabatic regime calculated within the small-q approach 
(solid lines) as function of Qc. In an exact theory the 
solid lines would end at Qc = 1 in the filled circles, rep- 
resenting the nonadiabatic local theory. The discrepancy 
is shown to be quite small and almost independent of the 
nonadiabatic parameter. In general, we can identify two 
different regimes: the first one, for small Qc's, is relevant 
for strong correlated systems where the critical tempera- 
ture is effectively enhanced with respect the conventional 
theory by nonadiabatic effects; the second one which rep- 
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FIG. 2. Tc (panel a) and Kfr (panel b) calculated by the 
small-q expansion of the vertex corrections for Qc = 1- Solid 
and dashed lines as defined in the previous caption. 

resent weakly correlated compounds where nonadiabatic 
vertex corrections produces a significant reduction of T^. 
While the first regime is expected to be well described by 
the small-q approach, Fig. ^ shows that is works qualita- 
tively well also in uncorrelated systems providing there- 
fore a unique tool to evaluate the nonadiabatic effects 
on the critical temperature. This appears even more im- 
portant since the analytical study of the vertex function 
in small-q expansion allows to understand from a mi- 
croscopic point of view the way through the which the 
nonadiabatic vertex function enhances or reduces Tc re- 
spectively in correlated and uncorrelated systems. 

In conclusion, we have investigated the nonadiabatic 
effects in uncorrelated systems. We have shown that the 
breakdown of Migdal's theorem and the consequent in- 
clusion of vertex diagrams lead to a strong suppression of 
the critical temperature Tc- Nonadiabatic effects can not 
be neglected in uncorrelated materials as well as in cor- 
related ones when the Fermi energy is comparable to the 
phonon frequencies. We show also that the nonadiabatic 
theory based on a small-q expansion, early introduced 
in previous papers, works quite well even in uncorreleted 
systems where no predominance of forward scattering is 
present. Of course, as already stressed above, this conclu- 
sion holds true as long as higher order vertex corrections 
can be neglected {XuipYi/Ep sufficiently small). 

As a final observation, we would like to remark that 
the present analysis has been restricted to the half-filling 
case. Away from half-filling, and in particular for Fermi 
levels very close to the bottom (top) of the band, the 
vertex function changes significantly its structure becom- 
ing roughly shapeless in momentum space and mainly 
positiveJ13 In such a situation nonadiabatic vertex cor- 
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FIG. 3. Tc as function of Qc as evaluated by the small-q 
expansion (solid lines). Dashed lines represents T in non 
crossing approximation and filled circles the nonadiabatic lo- 
cal theory for uncorrelated systems. From top to bottom: 
uo/Ef = 0.1,0.4,0.7 

rections can give rise|-tp an enhancement of Tc even for 
uncorrelated systems.til This case can be for instance rel- 
evant for the recentUz-. discovered superconductivity at 
Tc ~ 40 K in MgB2ji3 where the chemical potential is 
very close to the top of the tr-bands, with Ep ~ 0.5 eV, 
thus comparable to the phonon frequencies. 
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